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ABSTRACT 

Two types of models for the formation of semi-periodic concentric multiple 
shells (M-shells) around asymptotic giant branch (AGB) stars and in planetary 
nebulae are compared against observations. Models that attribute the M-shells 
to processes in an extended wind acceleration zone around AGB stars result in 
an optically thick acceleration zone, which reduces the acceleration efficiency 
in outer parts of the extended acceleration zone. This makes such models an 
unlikely explanation for the formation of M-shells. Models which attribute the 
M-shell to semi-periodic variation in one or more stellar properties are most 
compatible with observations. The only stellar variation models on time scales 
of ~ 50 — 1500 yrs that have been suggested are based on an assumed solar-like 
magnetic cycle. Although ad-hoc, the magnetic cycle assumption fits naturally 
into the increasingly popular view that magnetic activity plays a role in shaping 
the wind from upper AGB stars. 

Subject headings: circumstellar matter — planetary nebulae: general — stars 
AGB and post-AGB — stars: mass loss 

1. INTRODUCTION 

The arcs and rings which appear in the images of several planetary nebulae (PNe) and 
proto-PNe (Sahai et al. 1998b; Kwok, Su, & Hrivnak 1998; Su et al. 1998; Sahai et al. 
1999; Bond 2000; Hrivnak, Kwok, & Su 2001, hereafter HKS), as well as in one AGB star 
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(IRC+10216; Mauron & Huggins 1999, 2000), are thought to be concentric (more or less) 
semi-periodic shells; some shells are complete while other are not. Reviews of the arc and 
ring properties are given by HKS and Kwok, Su & Stoesz (2001; hereafter KSS). Not all 
observed arcs and rings are truly concentric, but in all cases their "center" is relatively (to 
their radius) close to the central star (KSS). I will use the fact that not all arcs are concentric 
to constrain models for their formation. Like Mauron & Huggins (2000) and to save space, 
I term these concentric semi-periodic shell multiple-shells, or M-shells. The time intervals 
between consecutive shells change from system to system and in some cases between shells 
in the same system, with typical time intervals of t s ~ 100 — 1000 yrs. 

Three classes of models were proposed for M-shells formation, (i) Binary models (e.g., 
Harpaz, Rappaport, & Soker 1997; Mastrodemos & Morris 1999) where the orbital period 
sets the time interval between the M-shells. The advantage of this kind of model is that the 
(semi)-periods are exactly in the orbital period range where the binary interaction, being 
tidal interaction (Harpaz et al. 1997) or accretion (Soker 2001), turns off and on along an 
eccentric orbit. However, there are strong arguments against binary models (e.g., Soker 
2000; Mauron & Huggins 2000). Hence, binary model will not be addressed in this paper, 
but I do allow for the possibility that a companion plays a passive role, e.g., spining-up 
the AGB star to initiate magnetic activity, (ii) Processes in the circumstellar matter, e.g., 
instabilities in dust-gas coupling (Deguchi 1997; Simis, Icke, & Dominik 2001). The relevant 
processes in this type of models occur in regions of size r a which are much larger than the 
stellar radius, r a 3> -R*. The basic advantage is that the dynamical time ~ r z J 2 / (GM) 1 / 2 , 
where M ~ 1M is the stellar mass, can be set to be of the order of the intershell time 
interval if r a ~ 20 — 100 AU. In the next section I argue that despite this advantage, there 
are problems with this class of models. (Hi) Semi-periodic variations in one or more of 
the stellar properties, leading to semi-periodic variation in the mass loss rate. Helium-shell 
flashes are ruled out (Sahai et al. 1998b; Kwok et al. 1998) because their typical inter-flash 
period is > 10 4 yrs, while pulsation periods are too short. This led me (Soker 2000) to 
suggest that a solar-like magnetic activity cycle in the progenitor AGB star is behind the 
formation of M-shells (see also Garcia-Segura, Lopez, & Franco 2001). 

In the present paper I compare models of the above types (ii) and (Hi) against properties 
of M-shells, and show that circumstellar models suffer from severe problems. The stellar semi- 
periodic variation can account for the M-shell properties, in particular if a binary companion 
is present. 
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2. CONSTRAINING AND SUPPORTING DIFFERENT MODELS 

2.1. Radiation Pressure and Opacity 

In this subsection I show that any model for the M-shell formation which is based 
on processes in the circumstellar matter is in conflict with wind acceleration by radiation 
pressure. The relevant process which forms the M-shells must occur before the wind reaches 
its high-Mach terminal speed. For simplicity I take a constant acceleration a over a radius 
r a much larger than the stellar radius r a R*. To form variation over a time interval t s , the 
acceleration time of each mass element should be of that order. The terminal wind velocity 
is Vt = at s , and the wind velocity in the acceleration zone R* < r < r a is v = (2ar) 1 / 2 = 
(2v t r /t s ) l l 2 . Assuming spherical symmetry, the wind density is 

_M_ Mil' 2 1 m 
9 4nr 2 v An{2v t yi 2 r^ K * < r - r a- UJ 

Since the material in the acceleration zone expands slowly, it is influenced by the gravity 
of the central star. The radiation pressure must overcome the gravity, since in the wind 
zone the pressure gradient is negligible. The total gravitational force on the material in the 
acceleration zone is 
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where equation (1) for p was used in the integration. Radiation can exert a force F ra ^ = (3L/c, 
where c is the speed of light, L the stellar luminosity and (3 is the effective average number 
of times a photon is scattered. Substituting typical values, e.g., for IRC+10216 (Winters, 
Dominik, & Sedlmayr 1994; Mauron & Huggins 2000) we find the constraints on the average 
number of times a photon is scattered in order to support the extended acceleration zone, 



Pa 




(3) 



In addition to supplying the counter force in the acceleration zone, the radiation supplies 
the momentum to the wind Mv t . The average number of scattering required here is 

^ = 0,74 (2 x 1O 4 L ) ( 2 x 1O~ 5 M yr- 1 ) (l5km s" 1 ) ' (4) 
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The total effective average number of times a photon is scattered is (3 t — Pa + Pw For these 
parameters f3 t — 3, which requires the wind to be optically thick. Indeed, the optical depth 
from a radius Rd < r < r a to infinity, where Rd — 10 AU is the dust formation radius, is 
given by 




where k is the dust opacity scaled with a typical value for AGB stars (Jura 1986; Winters et 
al. 2000), and the density was taken from equation (1). Only at r ~ 100 AU does the wind 
become optically thin. 

Simis et al. (2001) take L = 2.4 x 10 4 L Q and obtained average mass loss rate and 
velocity of M ~ 10~ 4 M Q yr _1 and v t ~ 40 km s -1 , respectively, which require f3 w = 8. 
This is a very high efficiency for wind acceleration, not found usually in AGB stars (Knapp 
1986). Such a large average number of times a photon is scattered requires that radiative 
transfer be followed through the optically thick acceleration zone. Simis et al. (2001) don't 
follow the radiative transfer, but rather assume optically thin envelope, making their results 
questionable. (From the density profile in the model of Simis et al. I find that there is a mass 
of ~ O.2M in the region between the stellar surface, R* = 6 AU, and the dust formation 
radius, Rd = 15 AU. This is more than the envelope mass in their stellar model. It is not 
clear what is the influence of this unreasonable assumption on their results.) 

To summarize, any model based on an extended acceleration zone requires high efficiency 
of radiation acceleration extending along the entire acceleration zone. However, for typical 
parameters, most, or even all, of the acceleration zone is optically thick, because the velocity 
is low and density high. Since radiation reaches outer region only after it is scattered to 
longer wavelength, acceleration tends to be less efficient in outer regions. I therefore consider 
unlikely that such models can explain the formation of M-shells. In any case, radiation 
transfer must be considered in building such models, as well as the stability of the process 
( Garcia- Segura et al. 2001). 



2.2. Non-concentric Arcs 

In three objects the M-shells are not concentric and/or intersecting each other (HKS): 
AFGL 2688 (the Egg Nebula), NGC 7027, and IRC+10216. Interestingly, the aspherical 
dense nebula in the center of each of these objects departs from axisymmetry: In the Egg 
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Nebula the eastern bright region ~ 6" from the center observed in the NICMOS image is 
more extended than the western bright region on the other side of the central star (Sahai et 
al. 1998a). The HST images of NGC 7027 (Latter et al. 2000) show a more extended and 
brighter region to the southwest side along the equatorial plane. The carbon star IRC+ 10216 
contains an inner elliptical structure on a scale of ~ 10" (Kastner & Weinntraub 1994). The 
variation in the spacing between the contours and polarization vectors (fig. 2 of Kastner & 
Waintraub 1994) are similar in nature to the variation between the M-shell spacing. The 
departure from axisymmetry appears also at ~ 0.1" from the center (Weigelt et al. 1998; 
Haniff & Buscher 1998). 

A plausible common explanation for both the non-concentric M-shells and the departure 
from axisymmetry of the inner bright region is a binary model (Soker & Rappaport 2001). 
The intersecting shells require the orbital period to be shorter, but not much shorter, than 
the intershell time interval. Each shell then, is ejected while the mass losing AGB star is 
at a different orbital phase, hence moving in a different direction. Since in all three objects 
the M-shells intersect close to the star, the orbital velocity of the AGB star around the 
center of mass v\ can't be too small (see Soker & Rappaport 2001). If two consecutive shells 
are ejected when the AGB star moves in opposite directions, i.e., half an orbital period 
apart, then the shell will collide when at a distance of ~ ARv t /(2vi) from the central 
star. In other cases they will collide further away. For these objects I find r\ < 5 — 10Ai?, 
hence v\ > 0.1 — 0.05t>t (v t is the terminal expansion velocity of the M-shells). I take as 
typical values v t = 15 km s -1 , and AGB and companion stellar masses of M 1 = O.8M , and 
M 2 = 0.6M Q respectively. The constraints on the orbital separation and period are (eqs. 
1-2 of Soker & Rappaport 2001) a < 100 AU and P or b ^ 10 3 yr, respectively. This binary 
model naturally explains the intersecting and/or non-concentric M-shells in models where 
the M-shells are formed during a short interval close to the AGB star. 

This binary model can't account for the non-concentric M-shells in extended acceleration 
zone models, where the formation time of each M-shell is about equal to the intershell time 
interval. The reason is that the required binary orbital period is of the order or shorter 
than the intershell time intervals. Hence, either each shell will be smeared, or else each shell 
have the average velocity, which is the binary center of mass velocity, and the shells will not 
intersect and not be displaced, i.e., they will be concentric. Therefore, another explanation 
should be found for the presence of non-concentric and/or intersecting M-shells in models 
based on an extended acceleration zone. 

I therefore predict that these three objects contain a binary companion of mass ~ 
0.5M Q , either a white dwarf or a main sequence star, with orbital periods inside the range 
of ~ 50 - 500 yrs. 



- 6- 



2.3. Polar-Disccontinueity 

In some systems there is a clear difference between the M-shell structure along the 
symmetry axis and that in the equatorial plane. This is usually attributed to different 
illuminating radiation, e.g., the searchlight beams in the Egg Nebula. By carefully examining 
the structure of some M-shells, I find the situation to be more complicated, (i) The different 
illuminating radiation can explain the brighter M-shell segments along the polar directions 
in the Egg Nebula. However, the M-shells in the Egg Nebula split, intersect, and possess 
irregular spacing along the polar direction much more so than closer to the equator. This 
is more prominent along the northern search light beams. This effect can't be attributed to 
different illumination, (ii) In IRAS 20028+3910 (HKS) there are three short arcs on each 
side along the symmetry axis, with the inner one being at ~ 1". However, there is faint 
nebulosity in other directions from the central star of IRAS 20028+3910, up to ~ 1.5" from 
the central star in the SW direction (Ueta, Meixner, & Bobrowsky 2000), where no arcs are 
observed (HKS). I therefore argue that the presence of arcs along the symmetry axis but 
not, or much less prominent, in other direction in IRAS 20028+3910, is real. 

The physical differences between equatorial and polar segments of M-shell that I argue 
for here are very tentative. In any case, it is interesting that magnetic activity, i.e., surface 
spots, are known to have polar and/or equatorial concentration, as well as uniform distribu- 
tion (e.g., Schrijver & Title 2001 and references therefin). Magnetic activity can naturally 
explain polar to equatorial variation (Gracia-Segura et al. 2011), if they turn out to be real. 

3. SUMMARY 

The goal of this paper is to examine two types of models for the formation of semi- 
periodic concentric multiple shells (termed here M-shells; also called arcs) against observa- 
tions. The comparison of models with observations suffers from the small number of objects 
with well defined M-shells. Binary models were not discussed here, as they were found in 
previous works to suffer from severe problems (Soker 2000; Mauron & Huggins 2000). In 
particular, the inter-shell time intervals are not constant. The binary model can be 'saved' 
only if another non-periodic component is added to the periodic orbital motion. I know of no 
mechanism to cause such a non-periodic variation component which can act together with 
the orbital period (beside the magnetic cycle which can act by itself; see below). 

In the present paper I find that models that attribute the M-shells to processes in an 
extended, up to ~ 100 AU, wind acceleration zone, result in an optically thick acceleration 
zone. This reduces the acceleration efficiency in the outer parts of the extended accelera- 
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tion zone, making such models unlikely. Simis et al. (2001) for example, assume that the 
acceleration zone is optically thin; hence, I consider their results questionable. 

Models which attribute the M-shell to semi-periodic variation in one or more stellar 
properties are most compatible with observations. The departure of the M-shells from being 
concentric in some systems can be easily explained by an orbital motion of the AGB star 
around the center of mass with a stellar companion. The stellar companion does not play a 
dynamical role in the formation process itself, but it may spin-up the AGB star to stimulate 
an efficient dynamo. The only stellar variation models on time scales of ~ 50 — 1500 yrs that 
have been suggested are based on a magnetic cycle (Soker 2000; Garcia-Segura 2001). From 
main sequence stars it is well known that in addition to a, more or less, uniform surface 
activity, magnetic fields, e.g., spots, can have polar and/or equatorial concentrations. This 
may explain the tentative claim I made here that in some objects there are structural differ- 
ences between M-shell segments along the polar directions and equatorial plane. Although 
the magnetic cycle is an ad-hoc assumption, it fits naturally into the increasingly popular 
view that magnetic activity plays a role in shaping the wind from upper AGB stars (e.g., 
Chevalier & Luo 1994; Rozyczka & Franco 1996; Soker 2000; Gardiner & Frank 2001). 

I thank Robert Link for comments on the original manuscript. This research was sup- 
ported in part by grants from the US-Israel Binational Science Foundation. 
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